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Abstract — The goal of this paper is to design compact support 
basis spline functions that best approximate a given filter (e.g., 
an ideal Lowpass filter). The optimum function is found by 
minimizing the least square problem (^2 norm of the difference 
between the desired and the approximated filters) by means 
of the calculus of variation; more precisely, the introduced 
splines give optimal filtering properties with respect to their 
time support interval. Both mathematical analysis and simulation 
results confirm the superiority of these splines. 

Index Terms — Spline, Interpolation, Filter Design 



I. Introduction 

THE conversion of continuous-time signals such as mul- 
timedia data with discrete and digitized samples is a 
common trend nowadays. This is mainly due to the existence 
of powerful tools in the discrete domain. However, the con- 
version of continuous-time signals into the discrete form by 
means of sampling may destroy all or some parts of the data. 
Under certain conditions on the continuous signal, such as 
bandlimitedness Q, the sampling process is guaranteed to 
be one to one; i.e., there should be a priori a continuous 
model. In spite of the technological movement toward digital 
signal processing, by the advances in wavelet theory Il2l-ll4ll. 
a revival of continuous-time modeling for the digital data has 
been triggered. Multiresolution analysis 0, (^i self- similarity 
171, (H, and singularity analysis 01 are inseparable from a 
continuous-time interpretation. It is therefore crucial to have 
efficient mathematical tools that allow easy switching from 
the digital domain to the continuous, and this is precisely the 
niche that splines, and, to some extent, wavelets, are trying to 
fill. 

In this field, polynomial splines, such as B-splines, are 
particularly popular, mainly due to their simplicity, compact 
support, and excellent approximation capabilities compared 
other methods. Spline-based methods have spread to various 
applications since the development of B-splines fTOl-fT?]. 

Though B-splines generate remarkable results in many 
applications, they are not the optimum solutions for filtering 
problems such as interpolation. This paper, focuses on the 
problem of designing optimal compact support splines which 
best approximate a given filter such as the ideal lowpass filter. 
In fact, the desired filter reflects the characteristics of the 
continuous-time model and can be arbitrary. 

The remainder of the paper is organized as follows: The 
next section briefly describes the spline interpolation method. 
In section [Till a novel scheme is proposed to produce new 
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optimized splines for interpolation regardless of the type of 
filtering. The performance of the proposed method is evaluated 
in section [IVlby comparing the interpolation results of the pro- 
posed method on standard test images to those of well-known 
interpolation techniques. Section |Vl concludes the paper. 

II. Preliminaries 
In this paper, the following notation and definitions are used: 

Definition 1. For a continuous-time signal x(t), a continuous- 
time signal Xp(t) and a discrete-time signal Xd[n\ are defined 
az follows, 

Xd[n]=<nT) (1) 



Xp{t) ^ x{t)p{t) = XdW{^-^T) 



(2) 



where p{t) = X^n^cx) ^{t — nT) is the periodic impulse train 
that is referred as the sampling function. ( Fig. [7]) 

The sampling period T = 1 is normalizes throughout the 
paper without any loss of generality. 

Definition 2. For a continuous-time signal x{t) and any odd 
integer m, x^{t) is a polynomial spline of order m if 

1) For any n G Z, x'^{t) would be a polynomial of the (at 
most) order m, in the interval [n, n + 1]. 

2) For any n e Z, x^{n) = Xd[n] (Interpolation property) 

3) x^ e C^-i(-(X),(X)) (Smoothness) 

According to the first property, m + 1th derivation of is 
zero in non-integer points, and is equal to an impulse train. 

Definition 3. For the polynomial spline x^ (t\ the polynomial 
spline coefficients x^ [n] are defined as, 



dt 



m+l ' 



\t) (3) 



To determine each polynomial of order m that is forming 
the x^(t), its m + l unknown coefficients should be found in 
order to satisfy the conditions 2 and 3 (Fig. O. If the goal is 
to discover a piecewise polynomial signal that is m — 1 times 
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Fig. 1. SampHng process modeled by multiplying an impulse train into a 
primary time continuous signal 
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differentiable with continuous derivatives, a natural way is to 
derive [n] according to Xd [n] and then calculate the integral 
of x'^{t), m + 1 times, i.e, 

xT{t) = I I x'^{to)dto...dtm-ldtm 
J — oo J — oo J — oo 

= (4) 
where u^{t) is the unity step function and for any /c G N, 

Proposition 1. If the ROC of X^{z) is not bounded by the 
unit circle (i.e, there exist z G i?OC{X^} such that \z\ > Ij, 
then x^{t) will be uniquly deriveble according to Xd[n]. And, 

xT{t)={{u"^+'*{u^+')-')*x,){t) (5) 

where {u^~^^)~^{t) is defined as the inverse ofu^^^{t), i.e, 
* ^^^^) W = d{t).AndXf{z) is the z-transform 

ofx^[n]. 
Proof: 

x^{t) = xT{t)p{t) 

= {u^+'*x^){t) 



(6) 



Hence, 



(7) 



The ROC of U^~^^{z) is \z\ > 1 and there is no zeros in 
this region either Since the ROC of X^{z) in not bounded 
by the unit circle, {U'^~^^)~^{z) and X^{z) have a region in 
common. Thus, 

x^(t)={{u^+^)-'*x^){t) (8) 
And according to (0), 

x^(t) = {u^^'^x^)(t) 

= (9) 



Definition 4. A discrete-time signal yd[n] is called an appro- 
priate signal if and only if it will be stable and have a unique 
and stable inverse yd^[n]. 



Sx (t)\n<t<n+i = a^[n]t^+a2[n]t^+ai[n]t^ +«o[«] 

SAn-) = SAn^) 

(n-) = (n^) 

dt dt 
d^S, d^S, 

(n-) = (n^) 

dt^ dt^ 



dS^ dS^ 

(n + r) = (n + l^) 

dt dt 

d^ d^ 

(„ + 1") = (n + 1+) 

dt^ dt^ 



Fig. 2. Spline of the order m conditions 
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Definition 5. For any continuous -time signal y{t), if yd[n] 
was an appropriate signal, then y{t) is defined as follows. 



m = {{yp)-'*y) it) 



(10) 



Proposition 2. y{t) is the impulse response of a filter with 
interpolation property, in the other word: 



yp{t) = s{t) 



(11) 



Proof: 



= yit)pit) 

= [{{yp)-'*y){t)]p{t) 

= {{yp)-'*yp){t) = 6{t) (12) 



Proposition 3. For any polynomial spline y^{t) ify^{t) was 
an appropriate signal, y'^if) is independent of y{t) and is 
only a function of m. 
Proof: 

yf{t) = {{{yT),)-'*yT)it) 
= {iyp)-'*y7)it) 

= (13) 
Hence y^(t) is only a function of m. ■ 

Definition 6. According to the above proposition c^{t) = 
y^(t) is defined as the cardinal spline of order m. 

From the equations ([5]) and ([T3l) it can be concluded that, 
the polynomial spline interpolation is a linear shift invarient 
process according to Xd[n] and can be exposed by c^{t). 
Where c^{t) itself can be deriven according to any arbitrary 
polynomial spline y^{t) that y^{t) is an appropriate signal, 
i.e. 



xTit) = (c-*x,)(t) 



(14) 



The above equation devides the whole interpolation process 
into a discrete-time and a continuous-time parts. If the y'^it) is 
choosen as a time limited basis, both parts of this process can 
be extremely simplified and a big amount of continuous-time 
calculation can be avoided. 

Proposition 4. Suppose that k is the least positive integer in 
which there exist a polynomial spline of order m like y^{t), 
that takes zero outside the interval (0, k) i.e. 



(15) 



then k = m -\- 1. 

Proof: Suppose that y^{t) satisfies the equation 4751 ) 
and Y^{z) is the z-transform of y^[n]. Where y^[n] is the 
polynomial spline coefficients signal ofy^{t) according to the 
definition It can be claimed that. 



{^-ir+i|i7(^-i) 



(16) 
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In order to prove f [7^ , define a sequence of polynomials 
{QnjJILo ^^^^ ^^^^ Qo{z) — ^(r(^~^) <^^dfor 1 < n < m, 

k 

Qn = z{^Qn-l) = J2 (17) 



n=0 



Also polynomial H is defined as follows, 



i=0 \n=0 / ^ ^ 

■ n=0 i=0 ^ ^ 

= -Ey™w(i-nr 



(18) 



n=0 



77zw5' according to the equation (0) /6>r t > k, H(t) is 
equal to y^{t), i.e, 



(19) 



Since H is a polynomial and is equals to zero for infinite 
amount of t, all of its coefficients are equal to zero, hence. 



= o^g™(i) = Q™_i(i) 



Qo(l)=0 (20) 



from the above equations it is concluded directly by induction 
that, 



(in . 



(21) 



Hence {z — On the other hand since 

= ELo^d^N^"' (ESI) cites that vXlm + 1] ^ 

hence, 

(m+1 — e,m+l+e) 

^ (22) 

thus k > m -\- 1. 

Finally it must be shown that there exist a polynomial spline 
of order m that is bounded by the interval (0, m + 1). Suppose 
that Y^{z) = {z-^ - 1)^+1 then. 



ynn] = {-iy 



m + 1 
n 

"m+l 



vTit) = u 



m+l 



Ey"N<^(i 



.n=0 



m+l 



(23) 
(24) 

(25) 



n=0 



Thus for all t > m -\- 1, y^it) = and the proof completed. 

■ 

Definition 7. The polynomial bspline of order m is defined 
as follows, 

m+l / ^ 1 \ 

r{t) ^ Y.^-lT[^^ \u'^^\t-n) (26) 

n=0 \ ^ / 



In order to have an FIR continuous-time calculation during 
polynomial spline interpolation process, it can be implemented 
by the polynomial bsplines, i.e, 



OO 

{t)= xMr{t-n) 



(27) 



(28) 



n= — OO 



III. Proposed Optimized B-Spline 

In many applications, it is desirable that the interpolation 
filter be depicted as an ideal filter, and the second and third 
conditions of definition [T] may not be important. In this section 
an optimized basis splines will be introduced to be replaced by 
the polynomial basis splines in order to have an interpolation 
process with the most possible coincidence with a desired 
filter. 

Definition 8. Let D denote the set of all continuous-time sig- 
nals that satisfy the dirichlet conditions, i.e for any y{t) G D 

1) y{t) have a finite number of extrema in any given 
interval. 

2) y{t) have a finite number of discontinuities in any given 
interval 

3) y{t) be absolutely integrable over a period. 

4) y{t) be bounded. 

First of all, an affine subspace of all signals that satisfy the 
dirichlet conditions will be defined, and then the optimized 
solution will be obtaind in this set by the calculus of variation. 

Definition 9. Let yd[Ti] be an appropriate signal that takes 
zero for all n < and n > m + 1, then X^iVd) is the set 
of all continuous -time signals y{t) that satisfy the following 
conditions, 

1) yeD 

2) Vn G N;y{n) = yd[n] 

3) Vt ^ (0,m + l);y(t) = 

Using y{t) G X^{yd) as a basis spline to interpolate Xd[n] 
according to the equations ([27]) and ([28b is a linear time 
invarient process with the impulse response y{t). 

Definition 10. The error function X^iVd) is defined 
as follows, 

/CO 
-co 



f 

J — c 

f 

J — C 



\H{{ypr'^y)^^p}-H^}\'df 



. (29) 

Where T is defined as the continuous time Fourier transform 
operator. 

Definition 11. According to the above definition, if be an 
appropriate signal that takes zero for all n < and n > m + l 
an optimized basis spline p^[x, p^] is defined as follows, 

p'^lx.p^]^ argmine,(y) (30) 
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Now, calculus of variation may be used in order to evaluate 
the optimum which minimizes the error ex{p'^). 

Proposition 5. Equation SSUi has a unique solution that 
satisfies the following property, 

* ^ * (p^)"' * (^)"'] * = * ^] * ^ (31) 

for all t G (0, m + 1). Where y(t) = y{-t) 

Proof: Considering 7 G X^(0)' variational derivation of 
ex{p^) ^ith respect to with j as a test function is equal 
to 



(e.(p^),7) 




dt (32) 



The proof of l[32\l is presented at box [1]. Since X^{p]i) 
boundless, in order to minimize ex{p'^), (e(p^), 7) should be 
zero for all 7 G X^(0)' which implies that the second term 
inside the integral should be zero for t G (0, m + 1), 



T{xp} 



Hp^} 



J^{Xp} - J^{x} 



(33) 



And this equation directly yealds 07]) . 



Thus, it is proven that the optimized basis spline which 
could give the best estimation of x, should satisfy (|3T1) . By 
defining 



v{t) ^ 



(34) 
(35) 



(I3TI) can be written as {v * p'^){t) — t(;(t)|^^(o,m+i)- Since 
this equation is only valid in a particular interval, cannot 
be used to obtain p^. But since v is an impulse train, p can 
be driven using matrix form, thus in order to derive from 
([3TI) , two sequences of functions are defined such that for any 
n G Z, 



Rnit) 
Wn{t) 




< t < 1 
o.w. 

< t < 1 
o.w. 

Now, (I3TI) could be written in matrix form as follows: 



^[0] 

^[1] 



-1] 



^[0] 
v[m\ v[m — l\ 



v[—m-\-l\ 
^[0] 





Rq 




r ^0 1 




Ri 




Wi 




- Rm - 




-Wm. 



(36) 



(37) 



(38) 



According to (l38l . {-Rn}™=i is derived and thus the optimized 
basis spUne is evaluated as 



n=0 



This optimized basis spline that finally calculated, preforms 
the least interpolating mean square error for x[n]. The basis 
spline could be calculated by the signal before sampling x, 




Fig. 3. Our Optimized Spline versus B-Spline, both of the order three. 
/3o{h,b} is the optimized basis spHne built for estimating ideal lowpass filter 
h{t) = £MM with b= (0.24,0.48,0.24) 



portion of x or statistic 's characteristics that are expected for 
X. Besides, smoothness of optimized basis spline and causality 
of prefilter are possible to be applied by setting /3s- 

Another application and even more important of equation 
(I3TI) is that by means of that it is possible to estimate any ideal 
interpolation filter by an optimized basis spline. In fact, spline 
dominance in against to other FIR windowed estimations is 
that an optimized basis spline with the same time interval can 
give more exact estimation for a desire ideal filter. 

Now the goal is to design p'^ such that p'^ would be the 
best estimation of h, which denotes the impulse response of 
a filter that has the interpolation property. 

Proposition 6. (Estimating an ideal filter) Considering h{t) 
as an impulse response satisfying the interpolation property 
and p^[n] as an appropriate signal, then 



argmin \\h - y\\^ = p^[h, p^] 



(40) 



Proof: 

eh{y) 



/oo 
\F{y*hp}-T{h}\^df 
-00 

/OO 
ITiyjJ'ihp} - T{h}fdf 
-00 

/OO 
\T{y} - T{h}\^df 
-00 



\\h-yh 



(41) 



Thus 



argmin \\h 



argmin eh{y) 



p^[h,p7] (42) 



Proposition 7. Consider h{t) as an impulse response with 
interpolation property and p^it) G X^{pT) ^ basis 
(39) spline which p^{t) best approximates h{t) over X^{pT) 

2- Then p^{t) satisfies the 



i.e, p^{t) = argmin II /i- 
following equation. 



(43) 



Proof: Follows directly from f l?71) and the fact that h{t) 
has interpolation property. ■ 

Fig. [3] shows the optimized B-SpHne built for estimating 
ideal lowpass filter h{t) = ^^^^ with pl{z) = 0.233z + 
0.480^2 + 0.233z^ and cubic B-Spline. And Fig. g] shows 
p^[h^p^]{t) in comparison to c^{t). 

IV. Simulation Results 

The performance of the proposed method for an ideal low- 
pass filter has been compared to the B-spline and the results 
are depicted in Figs. [3] and HI Fig. [3] shows the comparision of 
the optimized basis spline built for estimating an ideal lowpass 
filter and the cubic B-spline. Fig. H] shows L/^^ as compared to 
Ls. The optimized spline is superior to the B-spline method. 
The SNR values of these methods are 20.39dB and 13.15dB 
for the proposed method and the B-spline method, respectively, 
for m = 3. 

To consider practical applications, the method was tested on 
several standard monochrome images. These images are down- 
sampled to provide the low solution images for interpolation. 
In image applications, splines can be used for zooming and 
enlargements. For comparison, three other image interpolation 
methods are also simulated: 1 -bicubic interpolation, 2-wavelet- 
domain zero padding cycle- spinning [iT] and 3 -soft-decision 
estimation technique for adaptive image interpolation |[T5l . 
Table HV] shows the Peak Signal-to-Noise Ratio (PSNR) perfor- 
mance of these three methods when applied to the seven well- 
known test images. In all cases, the proposed optimized spline 
interpolation algorithm performed best among all methods. For 
high frequency content images, such as Barbara and Baboon, 
the proposed algorithm outperforms other methods by IdB. 

Since PSNR is an average quality measure, the spatial 
locations where the proposed algorithm produces significantly 
smaller interpolation errors than the other competing methods 
are plotted in Fig. [Jl The differences are more noticeable 
around the edge of the hat. The result of the present study 
compare favorably both subjectively and objectively. In addi- 
tion, a wavelet scheme based on cycle- spinning interpolation 
has been included to provide a comparison with a powerful 
method operating in the wavelet domain. 




-3-2-10 1 2 3 



Fig. 4. Comparison The performance of the ocomparisonur method and the 
cubic spline method for the ideal lowpass filter designing. 
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(e) (f) 



Fig. 5. Comparison of different methods for the Lena image: (a) The original 
image, (b) bilinear interpolation, (c) bicubic Interpolation, (d) WZP Cycle- 
Spinning [141, (e) SAI [T5j, and (f) the proposed method. 



TABLE I 

PSNR (dB) Results of the Reconstructed Images by Various 
Methods (Image Enlargement from 256 x 256 to 512 x 512) 



Images 


Bicubic |13 | 


WZP-CS 1 14 1 


SAI |15 | 


Opt. Spline 


Lena 


30.13 


30.05 


30.88 


32.29 


Baboon 


21.34 


21.70 


22.09 


22.50 


Barbara 


23.32 


23.88 


23.71 


25.10 


Peppers 


28.61 


28.60 


28.91 


30.64 


Couple 


26.73 


26.86 


26.96 


27.91 


Bout 


26.93 


27.07 


27.63 


28.50 


Girl 


29.97 


30.20 


29.94 


30.90 



V. Conclusion 

This paper has introduced a method for optimizing a com- 
pact support interpolating spHne for approximating a given 
filter in the least square sense. In particular,it demonstrated a 
newly proposed method for approximating the ideal lowpass 
filter. The interpolation results obtained by this method are 
better than those obtained by the conventional solutions, such 
as B- splines. Simulation results show about IdB improvement 
in most of the cases. In the future, we plan to focus on 



the application of these optimized spHnes for non-uniform 
samphng for 1-D and 2-D signals. 
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